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Abstract 

The cross section for exclusive production of ir + ir~ and 7r°7r° meson pairs in ultrarelativistic 
heavy ion collisions is calculated for LHC energy ^snn = 3.5 TeV taking into account photon- 
photon mechanism. We concentrate on the production of large two-pion invariant masses where 
■ the mechanism of the elementary 77 — > tttt process is not fully understood. In order to include 

5^ ■ a size of nuclei we perform calculation in the impact-parameter equivalent photon approximation 

(EPA). Realistic charge densities are used to calculate charged form factor of 208 Pb nucleus and 
to generate photon fluxes associated with ultrarelativistic heavy ions. Sizeable cross sections are 
obtained that can be measured at LHC. The cross section for elementary 77 — > tttt is calculated 
in the framework of pQCD Brodsky-Lepage (BL) mechanism with the distribution amplitude used 
to descibe recent data of the BABAR collaboration on pion transition form factor, using hand- 
bag mechanism advocated to describe recent Belle data as well as t and u-channel meson/reggeon 
exchanges. We present distributions in two-pion invariant mass as well as the pion pair rapidity 
Ph I for the nuclear process. 
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I. INTRODUCTION 



It was shown in several review articles [1] that the ultrarelativistic collisions of heavy ions 
provide a nice opportunity to study photon-photon collisions. This is due to the enhance- 
ment caused by the large charge of the colliding ions. Parametrically the cross section is 
proportional to Z\Z\ which is a huge number. It was discussed recently that the inclusion 
of nuclei sizes as well as realistic charge distributions in nuclei lowers the cross section com- 
pared to the naive predictions. Recently we have studied the production of p°p° pairs of 
muonic pairs [3;], of heavy-quark heavy-antiquark [4j as well as DD meson pair [5j. 

In the present paper we wish to study probably the simplest to measure exclusive produc- 
tion of pionic pairs. The elementary processes 77 — > tt + ti~ and 77 — > 7r°7r° have been studied 
in detail in the past (see e.g. jof). While very low energies are the domain of the chiral per- 
turbation theory [7|, at the intermediate energies one has to include also pionic resonances in 



the s-channel as well t and w-channel exchanges [6j, |8H10| . At low dipion invariant masses a 
huge contribution could come from a competitive photon-pomeron (pomeron-photon) mech- 
anism of exclusive p° production and it subsequent decay. The cross section for this process 
is very large (see e.g. [11]). At even higher energies \/i > 2 GeV the mechanism of the 
reaction is not fully understood. Brodsky and Lepage made a first prediction of the leading- 



order pQCD [12J which was further studied e.g. in [13|, |14j]. In general, the predictions of 
the pQCD calculation lay below the experimental data measured at LEP (THj] and recently 
by the Belle collaboration [16]. The next-to-leading order calculation has been carried out 
only in Ref. [13] and their result is not able to describe the present experimental data. The 
pQCD amplitude for the 77 — > 7T7T reaction depends on the pion distribution amplitude. 
It was believed for already some time that the pion distribution amplitude is close to the 
asymptotic form (6 x(l — x)). This turned out to be inconsistent with recent results of the 
BABAR collaboration for the pion transition form factor -F 7 » 77r for large photon virtualities 
18 1 . The authors of Ref. [I9J used a new model of the distribution amplitude which can 
describe the BABAR data. We shall use this model for the 77 — > tttt reaction. 

Some time ago Diehl, Kroll and Vogt (DKV) suggested that a soft hand-bag mechanism 



may be the dominant mechanism [20[ for wide-angle scattering at intermediate energies. In 
this approach the normalization as well as energy dependence of the corresponding cross 
section are adjusted to the world-data on the 77 — > 7r + 7r~ production (20| . 

In the present paper first we show how the different mechanisms describe the elementary 
data. Next we present our predictions for the nucleus-nucleus collisions. We will show distri- 
butions in the dipion invariant mass as well as in the pion pair rapidity. These are quantities 
which can be easily calculated in the impact-parameter equivalent photon approximation 
(b-space EPA). 



II. ELEMENTARY CROSS SECTION FOR 77 vrvr 
A. Perturbative QCD approach 

Basic diagrams of the Brodsky and Lepage formalism are shown in Fig. [TJ The invariant 
amplitude for the initial helicities of two photons can be written as the following convolution: 



M (A 1; A 2 ) = [ q dx £ dy <p n (x, pi) T^ lAa (x, y, p 2 ) fa (y, pfj , (2.1) 
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FIG. 1: Feynman diagrams describing the 77 — > (qq)(qq) — > vr7r amplitude in the LO pQCD. 



where \i x = min (x, 1 — x) y / s(l — z 2 ), fi y = min (y, 1 — y) J s(l — z 2 ); z = cos 9 12]. We 
take the helicity dependent hard scattering amplitudes from Ref. [HI]. These scattering 
amplitudes are different for 7r + 7r~ and 7r°7r°. It was proposed in Ref. (HJ to exclude the 
region of small Mandelstam t and u variables by multiplying the pQCD amplitude (12. ip 
by an extra form factor which cuts off the soft regions which were taken into account in 
Ref. explicity by including meson exchanges. The following form of the form factor was 
proposed in (2lj : 
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where t m the maximal kinematically allowed values of t and u. A reg is a cut-off 

parameter expected to be of the order of 1 GeV. The distribution amplitudes are subjected 
to the ERBLpQCD evolution [H, 23]. The scale dependent quark distribution amplitude 
of the pion [2J, |25|] can be expanded in term of the Gegenbauer polynomials: 
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where the expansion coefficients (only even above) can be written as: 
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C^/ 2 denote the Gegenbauer polynomials, 



QCD 



r, = 3, Np is the number of active quarks and A is the QCD scale parameter. 



Different distribution amplitudes have been used in the past Urn. Wu and Huang 

proposed recently a new distribution amplitude (based on a certain light-cone wave 



function): 
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This pion distribution amplitude at the initial scale is controlled by the parameter B. It 
has been found that the BABAR data at low and high energy regions can be described by 
setting B to be around 0.6. This pion distribution amplitude is rather close to the well 
know Chernyak-Zhitnitsky [27| distribution amplitude ((p^cz — 30x(l — x)(2x — l) 2 ). In the 
following (Eq. E5J we shall use B = 0.6 and m q = 0.3 GeV. Then A = 16.62 GeV^ 1 and 
(3 = 0.745 GeV. f n above is the pion decay constant. 

The total (angle integrated) cross section for the process can be expressed in terms of the 
amplitude of the process discussed above as: 

cr (77 tttt) = / 4-64 2 ; w2 g \M {Xi, A 2 ) | 2 dz , (2.6) 
where the factor 4 is due to averaging over initial photon helicities. 



B. Hand-bag model 

The hand-bag model was proposed as an alternative for the leading term BL pQCD 
approach j20[. It is based on the philosophy that the present energies are not sufficient 
for the dominance of the leading pQCD terms. As in the case of BL pQCD the hand-bag 
approach applies at large Mandelstam variables s ~ — t ~ — u i. e. at large momentum 
transfers. Diehl, Kroll and Vogt presented a sketchy derivation [20[ obtaining that the 
angular dependence of the amplitude is oc 1/ sin 2 9. Then the cross section integrated over 
cos 9 from — cos 9 to cos 9 for a charged pion pairs takes the simple form: 

/ . + -\ 4vra 2 m /cos^o . 1, l + cos^V , , |2 f0 

a(77^7T + 7r J= — - T — + -ln- -\\R nn (s)\ . 2.7 

v 7 s \ sin On 2 1 — cos On I 



Additionally, the ratio of the cross section for the 7r°7r° process to the ir + 7T process doesn't 
depend on 9 and is \. The nonperturbative object R nn (s) describing transition from a quark 



pair to a meson pair cannot be calulated from first principles. In Ref. 20] the form factor 



was parametrized in terms of the valence and non- valence form factors as: 

Rnn (s) = ^-a u ( + i« s ( , (2.8) 



9s V s J 9s V s 
where the authors of (lo[ have chosen s = 9 GeV 2 . The a u , n u , a s and n s values found from 



the fit in Ref. [20J slightly depend on energy. For simplicity we have averaged these values 
and used in the present calculations: a u = 1.375 GeV 2 , n u = 0.4175, a s = 0.5025 GeV 2 and 
n s = 1.195. The hand-bag approach was criticised in Ref. (26| . 



C. Meson exchanges in t or u channels 

Since several mesons (p, u, ax, 02, 61) decay into 77r channels this means that t and/or 
u channel exchanges of their virtual (space-like) counterparts may be important for the 
77 — > 7r + 7T _ and 77 — > 7r°7r° reactions. As an example in the following we consider u 
exchange for the 7r°7r° channel, p meson exchange also contributes to this reaction but its 
contribution is much lower (the corresponding coupling constant is 3 times smaller than that 
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for u meson exchange and it enters here in the second power already in the amplitude). So 
far the exchange of the tensor mesons was not discussed in detail in the literature. 

The amplitude for the 77 — > tttt reaction via vector meson exchange can be calculated 
by means of standard Feynman rules assuming tensorial form of the W7 coupling. The 
corresponding coupling constant can be obtained by fitting V — > n^y decay width. A simple 
and compact formula for the omega meson exchange amplitude was presented e.g. in Ref. 
(Io| . It can be written as: 

a h 2 

M (A lf A 2 ) = (X t (Ax, A 2 ) + X u (X u A 2 )) , (2.9) 

16 

X t (Ai A 2 ) = ei ( Al ) e2 ( A2 ) i! ( s ~ u ) + m i} - 2s ( e i ( e 2 (A 2 )p 2 } F2 ^ ^ 

X u (A 1; A 2 ) = 61 (Al) 62 (A2) {U {s ~ t)+ mt} ~ 23 {ei (Al) Pl} {£2 (A2) P2} Fl (u) , (2.11) 

u — 

where the size of the radiative coupling was obtained from the ratiative decay u — > 7r°7. 



In contrast to Ref. [10(] we include also vertex form factors (F w (t), F w (u)) which take into 
account the extended nature of the particles involved off-shell effects as well as high-energy 
reggezation. Not including the form factors leads, in our opinion, to nonphysical results, 
especially at large energies. Above W > 1.5 GeV the so-calculated cross section would 
significantly exceed experimental data. This point was not disccused in the literature as 
previous analyses were limited to rather low energies where the problem was not visible 
(note a remark in Ref. 0]). 

Using a vector particle propagators at high energy is not sufficient and one has to include 
reggezation. This is included in our calculation by multiplying the t and/or u-exchange 
amplitudes by the extra energy dependent factors: 

F u (t/u) = exp f 1 ( - ) • (2-12) 




The 00 trajectory is parametrized as a(t/u) = 0.64 + 0.8 t/u [28| and s = 1 GeV is taken in 
further calculations. 



D. Results 

In Fig. |2] we show the predictions of the hand-bag approach (solid lines), reggeized u 
- exchange (dotted lines) and the Brodsky - Lepage pQCD approach (dashed lines) for 
angular distributions of the 77 — > 7r°7r° reaction for W = 2.02, 2.26, 3,05 3.95 GeV. The 
pQCD results have been calculated in the case when F^ CD = 1. The cut-off parameter A w 
in Eq. (I2.12p was taken to be = 1 GeV. The results of different calculation are confronted 
with the Belle data. For the energies of present experiments the pQCD result is well below 
the experimental data. As can be seen from the figure the u - exchange may play a role only 
at large | cos#|. The result of the hand-bag approach starts to describe the data at energies 
> 3 GeV. 

In Fig. [3] we compare the pQCD 77 — > mi cross section for the pion distribution ampli- 
tude with and without pQCD evolution. The effect of the pQCD evolution on the angle- 
integrated cross section is very small, practically negligible. The data correspond to limited 
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FIG. 2: Angular distributions for the 77 — > 7r u 7r u reaction. The experimental data are taken from 
Ref. 



angular ranges given in the figure. The data for the 77 — > tt + tt reaction are from the 



ALEPH [15(, Belle [16j, CELLO |30|, CLEO [31], Gamma |32j, Mark II |33| and VENUS 



34 1 Collaborations. For the 77 — > -k°tt° reaction we present the Belle 35j and Crystall Ball 
36[ data. 

In Fig. H]we show the predictions of the hand-bag approach [2(| together with modern 
experimental data. The predictions can be taken seriously above the resonance region, 



i.e. when 



2.5 GeV. The parameters of the hand-bag contribution were adjusted 



to somewhat older experimental data. One can see that the hand-bag approach, 
consistent with the 7r + 7r~ data, slightly overestimates the 7r°7r° data. 



while 



In Fig. [5] we show the ratio of the cross section for the 77 — > n 



°7r° 



process to that for 
result and the 
cos 9 the ratio 



20 



the 77 — > 7r + 7r~ process. The dashed line represents the hand-bag model 
solid lines is for the Brodsky-Lepage pQCD approach. For larger range of z 
is smaller which means that the ratio is z dependent. The ratio is practically independent 
of the collision energy. In the present calculations, the z-averaged ratio for | cos(9| < 0.6 is 
about 0.2. The experimental error bars for the ratio (only statistical) were obtained with 
the help of the following formula: 
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FIG. 3: The pQCD cross section for 77 — > -k + tt~ (left panel) and for 77 — > 7r°7r° (right panel) as a 
function of photon-photon energy. The solid lines show the results for evolved 4> w (x, fi 2 ) and the 
dashed lines are for (j> w (x,/j,q) where /Iq = 0.25 GeV 2 was chosen. 
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FIG. 4: The hand-bag contribution for 77 — > tt + it (left panel) and for 77 — > 7t 7t° (right panel) 
as a function of photon-photon energy. 

The experimental data points are in between the predictions of the BL pQCD approach and 
the hand-bag model which further clouds the situation. 



III. THE NUCLEAR CROSS SECTION FOR THE PION PAIR PRODUCTION 



In our opinion the equivalent photon approximation in the impact parameter space (b- 
space EPA) is the best suited approach for applications to the peripheral collisions of nuclei. 
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FIG. 5: Ratio of the cross section for the 77 — > 7r u 7r process to that for the 77 — > 7r + 7r~ process. 
The experimental data were obtained based on the original Belle Colaboration 16|, |35[ data as 
explained in the text. 




FIG. 6: The Feynman diagram illustrating the formation of the pion pair as a result of the 
peripheral nuclear collision. 

In this approach absorption effect can be taken into account easily by limiting impact pa- 
rameter b > R\ + R2 ~ 14 fm. This approach have been used recently in the calculation of 
the muon pairs or p°p° pairs. The details of the b-space EPA have been described in @, |3[. 
Below we present a useful and compact formula for calculating the total cross section for 
the considered process: 

a (PbPb ->• PbPbnn; W^) = J a (77 ->• vrvr; W 77 ) 9 (|bi - b 2 | - 2FL A ) 

xA^(w 1 ,b 1 )A^(w 2 ,b 2 )27rfed6d^d6 y ^ I dW / 77 dF , (3.1) 

where the quantities iV (u, b) can be interpreted as photon fluxes associated with each of the 
nucleus and b x ,b y are auxiliary quantities which have been introduced in The photon 
flux is expressed in terms of the charge form factor. 

In Fig. [7] we show the modulus of the charge form factor of the 208 Pb nucleus for realistic 
charge distribution. The oscillations are related to relatively sharp edge of the nucleus. 
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FIG. 7: The modulus of the charge form factor of the 208 Pb nucleus for realistic charge distribution. 

Let us come now to our predictions of the nuclear cross sections. In Fig. [8] we show 
distribution in the two-pion invariant mass which by the energy conservation is also the 
photon-photon subsystem energy. For this figure we have taken experimental limitations 
usually used for the 7T7t production in e + e~ collisions. In the same figure we show our results 
for the 77 collisions extracted from the e + e~ collisions together with the corresponding 
nuclear cross sections for tt + it~ (left panel) and 7r°7r° (right panel) production. We show 
the results for the standard BL pQCD approach and for the approach proposed in Ref. (2o[ 
where an extra form factor given by Eq. f)2.2p was used to remove nonperturbative regions 
of small-angle scattering described at low energy in terms of meson exchanges. One can see 
that a difference occurs only at small energies which is not the subject of the present analysis. 
Above a/saT/v > 3 GeV the two approaches coincide. By comparison of the elementary and 
nuclear cross sections we see a large enhancement of the order of 10 4 which is somewhat less 
than Z\Z\ one could expect from a naive counting. 

In the e + e~ collisions the cuts on z = cos# are usually different for tt + 7t~ than for 7r°7r°. 
In the left panel of Fig. M we show the nuclear cross section for the same cut on z. In 
the Brodsky-Lapage pQCD approach the cross section for 7t + tt~ production is about order 
of magnitude larger than that for the 7r°7r° production. This is very different than for the 
hand-bag approach where the ratio is just \. As already commented above one can trust 
the pQCD results only for not too small energies and not too small angles or equivalently 
for not too small transverse momenta of pions. In the right panel we compare results of the 
Brodsky-Lepage pQCD approach (solid line) and results of the hand-bag approach (dashed 
line). Here in order to ensure validity of the both approaches we have imposed extra cuts 
on pion transverse momenta (p t > 3 GeV). At lower energies (W < 14 GeV) the hand-bag 
cross section is bigger than the cross section for the Brodsky-Lepage pQCD for the 7r + 7r~ 
production and the situation reverses at higher enrgies. For the 7r°7r° production the hand- 
bag cross section is always bigger than the BL pQCD cross section in the shown energy 
range. In this case the measured cross sections are not too big but should be measurable. 

As shown before the hand-bag approach better describes the elementary cross section. 
Therefore the hand-bag approach is used to estimate nuclear cross section. In the left panel 
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FIG. 8: The nuclear (upper lines) and elementary (lower lines) cross section as a function of 
photon-photon subsystem energy W 77 in the b-space EPA within the BL pQCD approach for the 
elementary cross section with Wu-Huang distribution amplitude. The angular ranges in the figure 
caption correspond to experimental cuts. 




FIG. 9: The nuclear cross section as a function of the 77 subsystem energy for the PbPb — > 
PbPbir + ir~ (green lines) and for the PbPb — > PbPb-K°TT° (red lines) reactions calculated for | cos 6\ < 
0.8 (left panel) and with an extra cut-off on pion transverse momentum pt > 3 GeV (right panel). 

of Fig. [10] we show pion pair rapidity distributions for different cuts. We hope that this 
figure may be a useful estimate of the cross sections for possible future experiments. In the 
right panel of Fig. [10]we compare the results of the BL pQCD approach and of the hand-bag 
approach for p t > 3 GeV (which by kinematics is equivalent to W 77 > 6 GeV). This is a 
region which was not measured so far in the e + e~ collisions. Nuclear experiment in this 
region should therefore discriminate between the two approaches. One could measure either 
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FIG. 10: The pion pair rapidity distribution. The left panel shows the result for the hand-bag 
model for different kinematic regions and the right panel compares the results for BL pQCD and 
hand-bag approaches for p t > 3 GeV, i.e. region not accessible so far in the e + e~ collisions. 



integrated cross section with cuts as well as study the ratio for 7r°7r° to tt + 7v as a function 
of accessible kinematical variables. 



IV. CONCLUSIONS 

In the present paper we have discussed a possibility to study the 77 — > nn processes in 
ultraperipheral ultrarelativistic heavy-ion collisions. 

In the present paper we have concentrated on the large two-pion invariant masses. First, 
we show how different reaction mechansims describe the large photon-photon energy data. 
We have discussed the pQCD Brodsky-Lepage mechanism with the distribution amplitude 
used recently to describe the pion transition form factors measured by the BABAR collabo- 
ration. For comparison we have considered the soft hand-bag mechanism proposed by Diehl, 
Kroll and Vogt. In addition we have considered also t and u channel u meson exchanges. In 
our opinion the situation in the measured energy range ^/s 77 < 4 GeV is not clear. 

The elementary cross sections have been used to make predictions for the exclusive pro- 
duction of pionic pairs in heavy-ion collisions. In order to concentrate on the interesting 
region where the pQCD may apply we have imposed cuts on pion angles in the dipion 
center of mass and on the pion transverse momenta. In addition, this allows to get rid of 
the soft and resonance regions. In the present paper we have presented predictions for the 
present LHC energy ^s^n = 3.5 TeV. The distributions in the two-pion invariant mass and 
pion-pair rapidity have been calculated and shown. 

Both the STAR collaboration at RHIC and the ALICE collaboration at LHC could mea- 
sure the cross section for the exclusive tt + tt~ production not only in the perturbative region. 
The region of resonances can be measured already with low statistics. Since the cross sec- 
tion for large invariant masses is smaller it requires good statistics. Having the absolutely 
normalized cross sections is very important in this context. In general diffractive nuclear 



11 



photon-pomeron mechanism can also contribute to the discussed region. Such a process is 
naively enhanced in nuclear collisions only by the Z 2 factor compared to the Z 4 factor for 
the mechanism discussed here. A real comparison to future data will require inclusion of the 
mechanism too. This goes, however, beyond the scope of the present analysis and requires 
further development in understanding nuclear diffractive processes. This is on our list of the 
topics of interest. 
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